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Coastal erosion, as inferred by measuring bluff recession is correlated with local 
wave height at twelve sites along the Monterey Bay coastline. Bluff recession rates are 
established by applying precise photogrammetric techniques to a 44-year time series of 
aerial photographs. Wave heights are determined from the U:S. Army Corps of 
Engineers Wave [Information Studies spectral wave climatology, where deep water 
gravity waves are hindcast from historic wind fields at three-hour intervals from 
January 1956 to December 1975. The deep water spectra are refracted to shallow 
water spectra at a nominal depth of 4 m. An erosion model is developed for Monterey 
Bay where the erosion process is modelled as a non-linear function of the 4-m 


significant wave height 
R = [AHs* + BHs + C(Tide + 1.02Hs -Clifftoe)]/Beach Slope 


The coefficients A, B, C are computed from a least squares regression of the modelled 
and observed recession rate values. The erosion model provides a reasonable 
representation of the erosion process in Monterey Bay, where the mean standard error 
between observed and modeled erosion rates is £0.17 m/yr. Adjustment of the wave 
energy coefficient, A, allows tuning of the model for high and low wave energy 


locations. 


TABLE ORGENIENTS 


I. EROSION PROCESSES IN MONTEREY BAY ..........-e0eeeeeees 9 
A.  I[NTR@D UPTON ~ cane een ent 7 eee ag 9 
B. BACKGROUND Hoe 7eee....... nna 9 
C. THE NEED FOR COASTAL EROSION DATA .............-. 14 
D> <OBIEGRIVES -.tas eee ie ee 16 
1] PHOTOGRAMMETRY =). a oe 18 
A BNCKGROUND .... 00>... Rate re 18 
DSITE SHBECTIONE . eee. tei ieccwicnsearouaale 18 
C. POTENTIAL ERRORS ASSOCIATED WITH 
PHOTO GRIAIMviedRA meet | ene 19 
D. INSTRUMENT ACCURACY AND SET UP.............-000- 23 
E. MEASUREMENT PROCEDURE ..........e.cececcuceceuees 25 
F. FIELD SURVEY ase ee ee ae 28 
G. EROSION DATA PRESHNTATI@Riee. on ce. a ..cn see, 30 
II. WARWIES AND TIDES =. nace I a 51 
A. WAVE CLIMATOLOGY =.) 51 
B. WAVE REFRAGIIC\ ie en 53 
C. WAVE SET-UBIANIDERU N= Upon 20 ee 59 
D. TIDES ..¢200¢..s0c ee 61 
IV. AN EROSION MODEL FOR MONTEREY BAY.........0.0eeeeees 65 
A. MODEL PRESENTA DI@iyen nn ee 65 
B. DISCUSSIONIOR RESUijica=nan rn nee 67 
C. CONCLUSIONS 7-) a I ee 70 
APPENDIX A: PROGRAMUSISMING PHON@w vi aeEn se 73 
APPENDIX B: PROGRAM LISTENG#) (i cR® qa ee 75 


Peele, RECESSION MODEL RESULTS 220.00... ccc cease. 


Piet OF REFERENCES 


oo 8 8 @ @ @ @ © © ®@ @ © © © & © © © ee © ee eh ehh hmhUlUhhU PMmhUhMhUhlUhOhUh MhUCUhhUcFmhUhlUhhUhhUCUlhrhUhhUhhmhUhhUhhUhHh Hh hUh hUhlUhHhUhUh Hh 


UNO CDSS ORG BS LTE) S10] 9 RSPR en er ce 


Ol 
Il. 
IV. 


VI. 
VII. 


VIII. 


Pa 


Xl. 
All. 


LIST OF Vaiss 


SELECTED soi ale Vico Seyret. 5 6a. + ei 20 
EROSIONSRASE RESULTS FOR SAND CITY «..4.:.0. 0... 5 er s 
EROSION RATE RESULTS FORK SOGMRRROeRIPCINID «25... cee 38) 
EROSION RATE IRESULTS ROR. MEN MI). Soe ats oe. oe BY 
EROSION RATE RESULTS FORPRUN@OMBREACH (..2....2252..98 39 
EROSION RA Te RESUL 15 FO Ries eioen orienta tanec nl eae 4 
EROSION RATE RESULTS FOR NORV iesiess LANDING ......8 44 
EROSION RATE RESULTS FOR ZMUD@ wy Sein S) <2. -s ee 46 
ACADEMY fice gece ee ec. AM chee cia 49 
CALCULATED RECESS IG Nie Cie eatin 7 2 aa ee 68 
RECESSION MODEL RESULTS (1956-1975) oo e eee eee ee ee ul 
OBSERVED RECESSION SUtaWeY © i... cee 0 ee 72 


foe 
eZ 
Ne) 
1.4 
201 
ja 
on 
2.4 
ZS 
2.6 
Ze] 
ZS 
Zo 
2.10 
oe) | 
a 
2ebS 
2.14 
ols 
2.16 
7a Ma 
onl 
5) I 
See, 
3.4 
BG, 
on) 
4.1 


PSO iG Wi Ss 


Soule Saitaserzieounty (Griggs, 1985). ooo, cscs e en enue... 10 
I OUUMmnmomeerey COUNtY (Griggs, 1985) (24... . se. cece ee cea e ee ee 1] 
Weicramivilonterey COUNLY (OTIGGS, PI6]) wastes ccs c sete sees e eee 12 
Pittorietamsport im vionterey Bay Oradiwe, 1986)". ... Se... 13 
iqesionmkatesevaluation Sitessae......0s.... De OO, 6 Bde 22 
Fechcim Displacemiame("VlOllitt, JOSU) ewan + ose eee ss ee ewe dese ec eee om) 
Mere pieteemenmevioritt, 1980) f2 2. cee ects eet ee cee ee wes 24 
eet eter AUIOUY VOM! OU amines ces cess ss + oe She eas eee HOM 23 
Ae eS eNO OMOUOCOMPArAUOMMMn Ne jase 2a «see ee Sede ee ee ee ee 26 
PRAT ak ata UNS IsOEMNe IVIALK occ stangiee. «cae as SOs se ss 2 OE m5 
Magee casiMemment PeCWMIGUe 6.....45...006 4s sees see ee dee cee eens es 2 
Dame mnn mE MOSIOM aNGSUITS |... ceyni Gia > sles y oases e cs Glas ey cs ee ee om 34 
DOrmmimoLeO TOME LOSION RCSUIUS a.9w 65 es ie 64 eM eee sec ee ees 36 
Pe: LOKOTMINCSUIES .. 1. mM IE OSe a roe he se ek ee 38 
PEICOMMR ECE ETOSIOUN IMCSUIUS: 4 clin a5 cs be Caudle See eka ee a cece ee 40 
Pamiceevab Erosion Nesults (Umsmoothiedy. ...... 006s... eee ee ese ee 42 
Pp vomiiemeciovenosiomermesults((Smoothed)mmes Sos ss. ewe eee 43 
iINortheNiess Landing Enosion; Results ....2...6%..ceee.... ie ee 45 
Zimuaowski munes Erosion Results’ Unsmoothed) ...................5. 47 
Zmadiowsk Dumes Erosion Results (Smoothed) ............0. ccc e ee aee 48 
Neagiiterey Bay Academy Erosion Reswhts 2... 2.6.0.6. 08..-00000: eer 50 
Wears enase | Grid with Archive Vocations Indicated ................... 54 
Wy ioeehase 1] Grdwith Arcinve Locations Indicated ..............566: 56 
Beem cce (Om laude! — iOh05edZ, 42 e ee. ee es he ss 5 egw ee ea 58 
lar memace (Oli. Academy, £ = 0.09 117, ii. ee cede eee ee ete ee es 60 
Tis ET COLE NOV ONGE) deNs BE 0 Seer er 62 
‘is SUNS (SIE ACY Sh se RR EAR 8 8 2 64 
Peer TOmmvnOGe! rome <2. Ghee ceca se eee etek ede ee eee ee eee 69 


ACKNOWLEDGEMENTS 


I would hke to first thank Professors Edward B. Thornton and Stevens P. Tucker 
for having the trust to accept me as a thesis student. Their patience, creativity, and 
academic enthusiasm have made this thesis experience both enjoyable and meaningful. 

Furthermore, I would like to express my sincere graditude to: the State-of 
California Department of Boating and Waterways for funding this research; Ms. 
Donna Burych for spectral transformation of the WIS data, development of the 
numerical bathymetry for Monterey Bay, and numerous hours of programming 
assistance, Mr. Stan Stevens, Librarian of the University of California, Santa Cruz, for 
providing aerial photographs; CAPT Glen Schaefer for procurement of precise survey 
equipment, and guidance in all aspects of field surveying; Mr. Sven Forst for his 
assistance in the development and installation of electronic photogrammetric 
equipment; ENS Pedro Tsai and AG2 Marty Shy for their efforts in the field survey; 
and the people who vandalized my car and stole the first draft of this thesis allowing 
me to do an even better job the second time. 

This thesis is dedicated to my wife, Nancy, whose efforts as a professional 
photogrammetrist insured the quality of measured data, and whose unfailing love, 
support and encouragement provided a productive home environment in which to 


complete this thesis. 


Il. EROSION PROCESSES IN MONTEREY BAY 


A. INTRODUCTION 


Each year, the energy generated in vigorous Pacific storms is released upon the 
Monterey Bay coastline. As the massive, long period waves break in the surf zone, 
their power is transformed into the forces that drive littoral transport and coastal 
erosion. These two processes scour and transport hundreds of thousands of tons of 
sand annually within the Bay. The popularity of Monterey Bay has brought pressure 
to develop its coastline; yet, past developments in this energetic zone, without proper 


regard for the coastal dynamics, have resulted in grave economic consequences. 


B. BACKGROUND 


Much of Monterey Bay area, particularly south of the Pajaro River, consists of 
low lying coastal plains bordered by foothills on the east, and active beaches on the 
west. This geomorphology is graphically depicted for South Santa Cruz County and 
fh emn Wronterey Coungy in Figures 1.1, #22, and 1.3. Centuries of sedimentary 
discharge from the San Lorenzo, Pajaro and Salinas Rivers have provided the 
abundance of sand required to build the coastal plains and beaches. This discharge, 
combined with a prevailing northwesterly breeze, and generally southerly longshore 
transport, has generated expansive sand dunes along the inshore side of the beach. In 
recent geologic time, however, the trend is towards erosion and not accretion. Erosion 
occurs when powerful storm waves, coupled with high tides, are able to rush up the 
foreshore and scour away at the base of the sand dunes. The undermining of the cliff 
causes the face of the dunes to collapse and fall onto the berm. Successive waves 
entrain the slumped sediment and carry it out into the littoral transport zone where it 
is either deposited into a less energetic region, lost into the Monterey Canyon, 720,800 
m?/yr, Oradiwe (1986), or mined by one of Monterey’s three sand companies, which 
collectively mine on the order of 413,100 m?/yr. Sand transport in Monterey Bay 1s 


shown in Figure 1.4. 
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Figure 1.1 
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Figure 1.3 Central Monterey County (Griggs, 1985). 
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Figure 1.4 Littoral transport in Monterey Bay (Oradiwe, 1986). 
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C. THE NEED FOR COASTAL EROSION DATA 


The delightful climate and picturesque setting of Monterey and Santa Cruz 
Counties attract an ever expanding population. One manifestation of this population 
increase is mounting pressure to develop the coastal areas. Development of several 
areas along the Monterey Bay coastline without adequate knowledge of the coastal 
processes, has resulted in considerable monitary losses. 

Seacliff State Beach, at the northern extreme of our research area, is a well 
documented example of a location that experiences storm damage with cyclic 
regularity. The seawall at Seacliff State Beach was destroyed and rebuilt three times 
between 1940 and 1978. In 1980 a new seawall was constructed at a cost of 1.7 million 
dollars. The structure was engineered to last twenty years, but unfortunately, in the 
winter of 1983, storm waves and their associated flotsam battered the structure until 
one third of the seawall was completely destroyed, and heavy damage was inflicted on 
many nearby camping facilities. The estimated reconstruction costs were $740,000. 
(Griggs, 1985) 

Further south, the exclusive Pajaro Dunes condominium development was 
constructed in a region of active sand dunes. The development’s location is one that 
experiences large scale episodic erosion followed by gradual accretion. The erosion 
occurs when the tidal stage is in phase with southwesterly storm waves. Three of these 
episodic erosion events occured between 1969 and 1983. The erosion was severe 
enough to cut the foredune back to the foundations of the development. Thousands of 
tons of riprap were placed in front of the development as a temporary measure during 
the 1983 storm scason. The following summer, a permanent revetment fortress was 
built at a cost of several million dollars. (Griggs, 1985) 

The Moss Landing Marine Laboratory was constructed on a site that has some 
interannual variability; however, accretion is the dominant morphologic process 
occuring on longer time scales. The foundation for the laboratory was laid in 
mid-December, 1983; yet less than one month later, storm waves cut back the foredune 
and encroached upon the construction site. Again, tons of riprap were implanted as a 
temporary measure, and a permanent seawall was installed the subsequent summer. 

Perhaps the most dramatic and well documented case of episodic erosion has 
occurred at Stillwell Hall on the Fort Ord Military Reservation. Studies by Lima and 
Sklavidis (1984) found the average cliff recession rate to be 1.83 m/yr. This dramatic 
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erosion is consistent with the refractive focusing of wave energy to the area (Chapter 
III). The U.S. Army has recently spent $400,000 dollars as a temporary measure, while 
it contemplates whether to spend millions for permanent protection, move this soldier’s 
recreation hall to a different site , or build a new one elsewhere. 

Lima and Sklavidis (1984) used precise, non-automated photogrammetric 
techniques to establish erosion rates at four locations in southern Monterey Bay: NPS 
Beach Lab, Monterey Sand Company, Phillips Petroleum property, and Stillwell Hall. 
Their results for the Phillips Petroleum property strongly corroborated the results of 
Thompson (1981), and their method of inferring coastal erosion through bluff recession 
(Thornton et. al., 1985) minimized the variance in photogranimetric measurements. 

Thompson (1981) viewed aerial photographs with a mirror steroscope to observe 
cliff toe recession at the Phillips Petroleum property in Monterey from 1939 to 1978. 
He concluded that the average recession rate was 0.6 m/yr. Moffitt (1968) used aerial 
photography to establish erosion rates in the vicinity of the Monterey Sand Company. 
Unfortunately, Moffitt referenced his data from the water line, the extreme variability 
of which resulted in rather large variance values. Nevertheless, he documented an 
average erosion rate of 2.1 m/yr. 

Lima and Sklavidis (1984) found cliff recession rates of 1.5 to 2.7 m/yr along a 
mile stretch of Sand City beach on either side of the Monterey Sand Company. They 
attribute this locally high erosion rate to the sand mining operations. 

The Monterey Beach Inn in Seaside and the Ocean House apartment complex in 
Monterey have both recently been threatened by major erosion events. The Monterey 
Beach Inn, previously named the Monterey Holiday Inn, constructed a permanent 
seawall to protect its investment, which at least thus far has been adequate. However, 
the Ocean House apartment management dumped five thousand tons of granite 
boulders onto the City beachfront adjacent to their complex as an emergency safegard 
measure. They are now required to remove the revetment from the public beach. 

The recent erosion events in Southern Monterey Bay may be natural phenomena 
in the annual and interannual variability of the coastline; or, as has been suggested by 
Griggs (1985), the increased erosion may be a result of sand mining in the Southern 
Bay. The reduction of sand availability caused by mining operations may inhibit the 
ability of large sand bars to form offshore during the winter. The annual formation of 
the offshore bars acts as a natural buffer against the powerful winter storm waves. The 


big waves break on these offshore sand bars and a majority of their energy is dissipated 


within the surf zone. Without the bars, or if they are reduced in size, energetic winter 
waves can expend their force directly on the beachface, and provide the energy required 
for a major erosion event. | 

The fact that several of the developed areas along Monterey Bay’s coastline have 
suffered the severe consequences of building in the dynamic coastal zone should be 
taken as a forewarning. The increasing population of the Bay Area will bring more 
pressure to further develop this delicate area. It is clear that quantifiable erosion rate 
data and wave climatology must play a key role in site selection and structure 
engineering in order to prevent the inevitable destruction of improperly situated, or 


engineered sites. 


D. OBJECTIVES 


The objective of this thesis is to develop a predictive model for coastal erosion in 
Monterey Bay. Erosion rates are inferred from bluff recession, which is measured using 
precise digital photogrammetric techniques from a 44-yr time series of aerial 
photographs. ‘Reference points that define the beach profile must be chosen to insure 
both sharp definition of the reference on an aerial photograph, and accurate depiction 
of permanent, long-term erosion. The variability of the foreshore profile caused by 
seasonal sediment transport renders it unsuitable as a measurement reference. 
Similarily, cliff material that has fallen onto the cliff toe makes that potential reference 
difficult to define, and generally unreliable. The sharp stereo presentation of the cliff 
top however, offers an ideal profile, and one that is not subject to short time-scale 
variability. Consequently, changes in the cliff top profile reflect permanent recession. 
Twelve sites have been chosen from Monterey, north to La Selva Beach and bluff 
recession is measured along 400-to 1300-m coastal intervals centered on each of the 
twelve sites. The resultant erosion rate data is used along with a 20-yr spectral wave 
climatology and tidal information in the development of a coastal erosion model. 

A majority of the previous erosion studies that led to quantifiable results used 
imprecise methodology in photographic measurement. These inadequate procedures 
produced questionable results and unquantifiable errors. The precise photogrammetric 
technique used by Lima and Sklavidis (1984) produced results with acceptable 
confidence limits. The manual mode of the stereocomparator that they used, however, 
induced operator measurement errors (reading the vernier), and was extremely time 


consuming. Hence only four sites were completed. 
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Our research employs digital equipment and techniques to minimize and make it 
easier to quantify operator and position measurement error. The variance in bluff 
recession data is minimized through numerical integration of the chosen coastal 
segments, and a comprehensive survey of the twelve coastal sites insures that the model 


will be representative of the erosion processes in Monterey Bay. 


Vi 


Il. PHOTOGRAMMETRY 


A. BACKGROUND 


Precise photogrammetry is the currently recognized and practiced cartographic 
method of producing accurate, cost effective, topographic maps and _ shoreline 
manuscripts for nautical charts. Its accuracy for measuring events occuring on long 
time scales, such as cliff recession, is far superior to the previously employed methads 
of inferring coastal erosion through comparison of historical maps and charts, or 
subjective “hand” measurements on aerial photographs. 

This research uses planimetric and precise metric photogrammetric techniques to 
establish cliff recession rates. Planimetric photogrammetry is employed to qualitatively 
survey the selected measurement sites. At this initial stage, the photogrammetrist 
views the photo model stereoscopically to select the scaling points, identify local photo 
control points, determine the clarity of the cliff line, and establish the overall 
servicability of the photographs. Once these points have been identified, a field survey 
is performed to obtain the ground distance between scaling points, the height of the 
berm, and the elevation of the cliff top (above MLLW). Accurate quantitative 


measurements are then made using a precise stereocomparator. 


B. SITE SELECTION 


In any erosion study, it is advantageous to establish the erosion rates for the 
entire coastal survey area; unfortunately, the use of such a method on large scale 
erosion problems becomes economically unfeasible. A viable alternative, however, is to 
choose a number of study sites in the survey area that are representative of the overall 
erosion processes. This thesis examines twelve sites in Monterey Bay from the City of 
Monterey north to La Selva beach. Erosion information has been calculated for eight 
of the twelve sites examined, and results from Lima and Sklavidis (1984) are used for 


the remaining four sites. The specific criteria for site selection are listed below: 


* Is the site of economic interest ? 


* Is the site of scientific intcrest ? 
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* Is the site of political interest ? 


* Does the site lend itself to the use of photogrammetric techniques ? 


* Will the site provide representative erosion rates ? 


The total number of sites had to be limited to what was felt could be accomplished 
within the given time constraints. Additionally, the sites had to be spaced along the 
Bay's coastline in such a manner as to provide representative coverage, and all sites 
had to be photogrammetrically usable. The geographic positions of each of the sites 
and their associated names are provided in table I. The actual study area consisted of 
a 500-to 1300-m stretch of bluff in front of each site. The site locations are shown in 
migube 2.1. 


C. POTENTIAL ERRORS ASSOCIATED WITH PHOTOGRAMMETRY 


The primary sources of error in photogrammetric data are caused by image 
displacement and scale variation. Specifically, photogrammetric errors can be classified 
into three basic catagories: relief displacement, tilt displacement, or terrain induced 
scale variations. 

Relief displacement in a vertical photograph occurs when an image is displaced 
from its actual position due to an elevation difference between the image and the 
photographic principal point. Images of points on photographs with significant 
topographic relief will be displaced radially outward and images at lower elevations 
radially inward from the photographic principal point, as shown in figure 2.2. 
Displacement has the effect of portraying distances between images at higher elevations 
to be greater than they actually are, and distances between images at lower elevations 


to be less than they actually are. Relief displacement can be quantified and 1s given by 
Com (F< hh) eR) 


where d is the radial displacement distance, r is the radial distance from the principal 

point to the image, h is the image elevation above datum, and H is the flying height. 
Tilt displacement occurs when an image is displaced from its actual position 

resulting from a tilt in the vertical camera angle at the time of photo exposure. Ona 


vertical photograph, the nadir, n, and the principal point, 0, are indeed the same point. 
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Station Name 


NPS Beach Lab 


Philips Petroleum 


Sand City 


Monterey Sand Co. 


South Fort Ord 


Stillwell Hall 


Marina 


Rincon Beach 


Marine Lab 


North Moss 


Zmudowski Dunes 


Monterey Bay Academy 


TABLE I 
SELECTED SURRY Ey shies 


Latitude 


36° 36'10" 
36° 36°36" 
36° 37°00" 
Someta 0) 
36° 38'05” 
36° 39°40" 
Bormalvoss 
36° 46'30” 
36 °47°35” 
36°48'25” 
36°49°40” 


36° 54’26” 
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Longitude 
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Iies0 5a: 


1D e40nss: 


121°49'17” 


er 48°32” 


47 43” 


Peay 32" 


LQeAT2Ss 


Le 507 


Ieis0 55c 


W 


W 


W 


W 


W 


W 


W 


W 


However, when a camera is tilted at the time of exposure, its bore sight will describe 
the photograph’s principal point, and a plumb line from the cainera’s optical center 
will describe its nadir, obviously now, not the same point. This distinction is important 
because image displacement on a tilted photograph will occur along radials of the nadir 
point, and not the principal point as was the case with relief displacement. 
Photographs with tilt displacement can be rectified by projecting their coordinates onto 
a new plane that is parallel to the map plane, or datum. The coordinates projected 
onto the new plane then comprise the equivalent vertical photograph. Tilt information 
was available for all of the photographs selected for this research, and only 
photographs with less than three degrees of tilt were used to establish cliff recession 
rates. The image displacement resulting from a 3° tilt of the horizontal flight line axis 
is less than the photocomparator’s resolution ability; thus, all photographs employed in 
this thesis were considered vertical. The geometry of tilt displacement is given in 
Picuue 26, 

Terrain-induced scale variation is another effect of vertical relief on a 
two-dimensional photograph. A manifestation of scale variation is that higher 
elevations will appear on the photograph at larger’scales and lower elevations will 
appear at smaller scales (assuming a constant flying height). This variation is shown in 
Figure 2.4. In a photogrammetric mission, scale variation can be decreased by 
increasing the flying height and camera focal length. 

There are two procedures for minimizing the effect of the aforementioned errors 
in a photogrammetric data base. The first procedure is to apply point specific tilt, 
relief, and scale correctors to each measurement taken on the photograph. 
Unfortunately, hundreds of measurements are typically taken on each photograph, and 
this process is slow and laborious. The second procedure, and the method adopted by 
this research, is to choose ground control, and scaling points that are both, as close as 
possible to the measurement region, and at nearly the same elevation. Proper choice of 
the control and scaling points insures that the photogrammetric measurements are 
made over distances close to the baseline, and that the photoscale is highly accurate for 


the measurement region. Photoscale is given by 


S = (ab)/(AB) (2.2) 
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Figure 2.2 Relief Displacement (Moffitt, 1980). 


‘where AB is the measured ground distance from the field survey, and ab is the 
photodistance measured with a photocomparator. A more detailed discussion of 


photogrammetric errors and their correctors can be found in Moffitt (1980). 


D. INSTRUMENT ACCURACY AND SET UP 


Photogrammetric measurements were made on a modified Zeiss stereocomparator 
interfaced to an IBM Personai Computer (IBM-PC). A schematic of the Zeiss 
instrument is provided in Figure 2.5. The digital modification was accomplished by 
placing Teledyne-Gurley DMC-400 encoders on the stereocomparator’s X, Y, and 
X-parallax adjustment controls. Each encoder outputs 8000 counts per one revolution 
of the comparator’s controlling screws. Horizontal coordinates can be read from the 
instrument’s verniers to an accuracy of 0.01 mm, and the X - parallax control can be 
read to 0.001 mm (Hallert, 1960). For calibration, numerous readings were taken 
along each of the instrument’s three verniers, and paired with their associated encoder 


counts. A linear regression was then performed for each of the three paired data sets 
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Figure 2.3 Tilt Displacement (Moffitt, 1980). 


to provide encoder count to horizontal distance transfer functions. These functions 
were evaluated for uncertainty by first calculating the standard deviation, or standard 
error for each axis, which were found to be: De 0.009 mm, Cy = 0.024 mm. The 
RMS radius of the error circle, or position error, a , for an uncorrelated, orthogonal 
coordinate system is given by Heinzen (1980) as 
d= (67, + 6° )!/2 (2.3) 

and 1s found to be 0.022 mm. This photo position error can be converted to ground 
position error by multiplying d. by the photoscale. For example, the ground position 
error on a 1:12,000 scale photograph would be two tenths of a meter. 

Stereo vision presents the sharpest possible presentation of the standard 9” x 9” 


photographic data. A pair of overlapping photographs are placed on plates 1 and 2 of 
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Figure 2.4 Scale Variation (Moffitt, 1980). 


the photocomparator (Fig. 2.5), each of which may be rotated in a plane normal to the 
vertical axis of the photograph. The image on plate 1 is rotated until the X-fiducials 
are parallel to the flight line (X-axis of the photocomparator). The image on plate 2 is 
also rotated to align the X-fiducials with the flight path, but additionally, a Y-parallax 
adjustment must be made to bring the flight line on image 2 into coincidence with the 
flight line on image 1. The Y-parallax will have been completely removed once the two 
flight lines are in coincidence, and the operator will be able to see two marks in the 
stereo viewer. One mark, on image 1, will be stationary, the other mark will appear to 
float above the photograph ( Figure 2.6). The X-parallax or elevation, control is then 
adjusted until the floating mark becomes coincident with the fixed mark. The 
conjugate points on the two photagraphs are now viewed as exactly the same point 


and precise stereo measurements can be made. 


E. MEASUREMENT PROCEDURE 


Several points are selected that appear on each stero model in a particular 
photographic time series. These points are: a cocrdinate system origin, O; a baseline 
definition point, B; and numerical integration limits, L; and L,, (Figure 2.7). Annual 


and interannual variability of the sand dunes render it impossible to repeatedly locate 
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Figure 2.5 Zeiss Sterophotocomparator. 


the integration limits by sight on a photograph; instead, these points ‘are accurately 
located by the intersection of the cliff edge and a line normal to the cliff that has been 
defined by two reliable points, IP, and IP, discernable on the photographs. The same 
scaling points are typically selected for each stereo model. This procedure is not 
always possible, however, and additional scaling points sometimes have to be chosen. 
Such additional points must be supported by supplimental field work. 

Digital coordinates are obtained through employment of an interactive program 
which is designed to communicate between the comparator’s encoders and the 
IBM-PC. The program, PHOTODAT, accepts commands from the IBM-PC keyboard 
and writes the desired coordinates into a station data file. PHOTODAT is listed in 
Appendix A. The digital data acquisition procedure is very specific and requires the 
origin of the photographic coordinate system be defined first. This step is 
accomplished by moving the comparator mark to the designated position, O, and 
entering it. The operator must then move the mark to the baseline definition point, B, 
and enter its position. The X and Y values of point B define the coordinate system 


rotation angle @ given by 
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Figure 2.6 X-Parallax and the Floating Mark. 


@ = arctan(Y/X) (2.4) 


All points entered subsequent to the baseline definition point are projected onto a new 


coordinate system via the coordinate system rotation matrix 


i [ cosa ine | (2.5) 
= A 
¥ —_ cosa ne 


The scaling parameters are the next points entered, and the first points to be projected 
onto the new coordinate system. Once the photographic scale has been accurately 
defined, bluff measurenents can begin. The cliff edge is outlined by a dense network of 
points entered monotonically along the chff from point L, to L,. The actual cliff 
outline is approximated by line segments whose lengths vary slightly with topography 
but are typically less than 0.25 mm at photoscale. A trapazoidal rule integration was 


performed to determine the area between the cliff edge and the baseline using 
: 
DE CY; + Yi4p/2 * (Xi41 - XG] (2.6) 
url 


wy 


Areas for adjacent years in the time series were subtracted to provide the actual eroded 
area, and that figure was multiplied by the mean cliff elevation to determine the eroded 
volume. The actual coordinate system rotation, scaling, and numerical integration are 
performed by an interactive program, WRKPROG (Appendix B). WRKPROG 
retrieves raw data from the PHOTODAT files, accepts inputs from field surveys, and 


executes the numerical computations. 


F. FIELD SURVEY 


A separate field survey was performed for each of the eight erosion sites. One 
task of the field party was to determine the authenticity of the photo control and 
scaling points. If there was any doubt as to whether a particular reference point had 
been moved or changed sometime during the 44-year photographic time series, that 
point was rejected, and another selected. Additionally, accurate distances had to be 
measured between the ground scaling points to determine the precise photographic 
scale. This operation was accomplished using a Keuffel and Esser, Ranger-IV, laser 
distance measuring device. The K&E Ranger-IV has an accuracy of £5 mm + 2 ppm 
(Keuffel and Esser, 1975). At typical measurement distances of 800 m, this uncertainty 
translates to £5 X 10° mmona 1:12,000 scale photograph. A standard deviation of 
+5 x 10% mm is one full order of magnitude less than the Zeiss photocomparator’s 
resolving power; consequently, no accuracy was sacrificed through use of the 
Ranger-I1V. The scaling distances at each station were measured ten times and 
averaged to insure the best possible results. Temperature, pressure and humidity were 
observed at each measurement site, and the measured distance was corrected for 
atmospheric refraction effects. 

The elevation of the cliff toe was measured at 75-m intervals along each site’s 
beachfront. The elevation measurement was observed using a surveyor’s level and a 
Philadelphia rod. The actual cliff toe elevation was leveled from the mean water line, 
and the measurements were subsequently reduced to datum (MLLW) using tidal 


correctors. 
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Figure 2.7 Photo Measurement Technique. 
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G. EROSION DATA PRESENTATION 


Erosion results are presented for eight sites along Monterey Bay: Sand City, 
South Fort Ord, Marina, Rincon Beach, Moss Landing Marine Lab, North Moss 
Landing, Zmudowski Dunes, and Monterey Bay Academy. 

Typically, 75 to 150 points are digitized from each stereo model. Integrated areas 
for adjacent photo models in the time series are subtracted, and the difference is the 
recessed area for the intervening time period. The recessed area is divided by the 
longshore observation distance (length of the digitized bluff), and that value is divided 
by the number of years between the photo models. The result is a cliff recession rate 
in meters per year. Negative values in the tables below indicate recession, whereas 
positive values indicate accretion. A standard deviation is calculated for the recession 
rates at each site, and that value is stated in the respective site’s result table. The 


number of standard deviations by which a particular recession rate drifts away from the 


Cd a a 


mean value is listed in the tables’ “|x; - x|/o6 ” columnns, and provides a measure of 
temporal variability. Assuming episodic recession, a high value in the sigma column 
would signal an erosion event; however, it is unportant to also consider the length of 
the intervening time period. Long time intervals will tend to alias erosion events, and 
conversely, erosion events occuring within a short time period will be highlighted. 
Average erosion is inferred by multiplying the recession rate by the average cliff 
elevation. Not surprisingly, areas of refractive wave focusing resulting in energy 
convergence experienced the highest recession rates; whereas sites located at the head 
of Monterey Canyon, where there is refractive divergence, experienced long-term 
accretion. It is noteworthy that all sites experienced at least one significant erosion 
event between 1980 and 1984. 

Sand City is a site that does not have a well defined coastal bluff. The bluff was 
approximated by digitizing the crest of the seaward-most sand dune. While recession 1s 
the dominate process occurring over the 44-year time series, significant accretion 
occurred between 1940 and 1956. This point is interresting, because - whether this 
growth event was a manifestation of sand dune variability, or bonifide accretion - the 
process made an abrupt change late in the time period. In the years following 1966 
Sand City displayed recession rates comparable to those of Fort Ord, where high wave 
energy drives the highest recession rates recorded in Monterey Bay. The high recession 


rates observed at Sand City are not corroborated by the wave energy studies discussed 
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in Chapter III. The recession results at Sand City most likely reflect the presence of 
the sand mining industry which is most prevalent in that region of the Bay (Table II 
and Figure 2.8). 

The highest recession rates we observed were at South Fort Ord. The recession 
is depicted in Table III and Figure 2.9, where the average bluff receded fifty meters 
from 1940 to 1984. This high rate is consistent with the refractive wave energy 
focusing which appears to be a maximum in this region. 

Marina is also a site where high wave energy produces high recession rates. 
Additionally, numerous observations suggest that Marina is an area of high wind 
energy, perhaps due to geographic focusing of the seabreeze. The high wind energy 
and the variable nature of the non-solidified coastal bluff cause high variance in the 
observed recession rates, and may, in part, be responsible for the anomolous accretion 
observation in the 1970 time frame (Table IV and Figure 2.10). 

Station Rincon Beach at the northern “edge” of the South Bay’s energy 
convergence zone displays continuous recession from 1940 to 1984. The inference of a 
convergence zone boundary is made as wave energy increases south of Rincon Beach, 
and stations just north of Rincon exhibit long-term accretion. The recession at Rincon 
Beach is seen in Table V and Figure 2.11. 

Accretion at ‘station Marine Lab is evidenced by a shift in the entrance of 
Elkhorn Slough which occurred sometime between 1940 and 1966. This shift is clearly 
depicted in Figure 2.12. Figure 2.13 omits the slough entrance to present the 
recession/accretion events as they occurred after 1966. The numerical values are seen 
in Table VI. 

The low wave energy caused by refractive divergence at the head of Monterey 
Canyon is responsible for littoral deposition and the resultant accretion at station 
North Moss Landing (Table VII and Figure 2.14). 

Station Zmudowski Dunes appears to be the northern conjugate of Rincon 
Beach. It is a location where higher wave energy produces continuous recession, yet 
only 1.5 miles to the south, lower wave energy allows the accretion process to 
dominate at station North Moss Landing. Figure 2.15 graphically depicts the 
recession, and the recession rates are provided in Table VIII. The variability of 
McClusky Slough is evident in figure 2.15 and 1s omitted to present Figure 2.16. 

Recession is the dominate process occurring at Monterey Bay Academy between 


1956 and 1984. The wave energy focused to the north of Monterey Canyon is typically 
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not as high as energy which has been focused to the south. This point is almost 
Intuitive, as a majority of the winter swells originate from extratropical cyclones in the 
Gulf of Alaska. Accordingly, the recession rates at Monterey Bay Academy are less 
than those at its southern conjugate station South Fort Ord. The recession at 


Monterey Bay Academy 1s presented in Table IX and Figure 2.17. 
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TABLE I] 
EX@SION RATE RESULTS FOR SAND CITY 


STATION NAME: Sand City 

LATITUDE : 36°37'00" N 

LONGITUDE : 121°51’15” W 

MEAN CLIFF TOE ELEVATION (above MLLW) : 4.53 M. 
MEAN CLIFF TOP ELEVATION (above MLLW) : 10.04 M. 
LONGSHORE OBSERVATION DISTANCE : 540.00 M. 
BEACH SLOPE : 1:40.00 


SIGMA : 1.00 M. 


TIME INTERVAL RECESSION —[x;-X/ AVG. EROSION 
[M/YR] 0 ((Me/YR)Y/M] 
1940 - 1956 +0.47 1.44 + 4.72 
1956 - 1970 -2.30 1.33 -23.09 
1970 - 1976 a1.22 0.25 31225 
1980 - 1984 -1.34 0.37 Bees 
1940 - 1984 i -9,74 
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Figure 2.8 Sand City Erosion Results. 
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PASE Til 
EROS. Neate RESULTS FOR SOUTH FORT ORD 


STATION NAME: South Fort Ord 

LATITUDE : 36°38'05” N 

LONGITUDE : 121°49’55” W 

MEAN CLIFF TOE ELEVATION (above MLLW) : 6.19 M. 
MPA ClIFE TOP ELEVATION (aboveiMLLW) : 35.39 M. 
LONGSHORE OBSERVATION DISTANCE : 692.00 M. 
BEACH SLOPE : 1:36.36 


SIGMA : 1.47 M 





TIME INTERVAL RECESSION __x,-3 AVG. EROSION 
[M/YR] C ((M3/YR)/M] 

1940 - 1956 -1.09 0.02 -38.58 

1956 - 1966 0.01 0.75 0.35 

1966 - 1980 ahs 0.01 -39,99 

1980 - 1984 BG 1.86 -136.85 

1940 - 1984 (7 re -39.64 
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Figure 2.9 South Fort Ord Erosion Results. 
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WWevelble Ia 
EROSION RATE RESULTS FOR MARINA 


STATION NAME: Marina 

LATITUDE : 36°41'58” N 

LONGITUDE : 121° 48’32” W 

MEAN CLIFF TOE ELEVATION (above MLLW): 4.92 M. 
MEAN CLIFF TOP ELEVATION (above MLLW): 24.31 M. 
LONGSHORE OBSERVATION DISTANCE : 1283.00 M. 
BEACH SLOPE : 1:40.00 


SIGMA : 1.67 M. 





TIME INTERVAL RECESSION |x;-X| AVG. EROSION 
[M/YR] o [(M?/YR)/M] 

1940 - 1956 -0.99 0.36 -14,83 

1956 - 1970 -1.12 0.44 -27.23 

1970 - 1980 oO) 1.60 + Oo.0'l 

1980 - 1984 a We 1.05 -51.54 

1940 - 1984 SQGS 8 8  denade -9,24 
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Marina Erosion Results. 


Figure 2:10 
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TABLE V 
EROSION RATE RESULTS FOR RINCON BEACH 


STATION NAME: Rincon Beach 

LATITUDE : 36°46'30” N | 

LONGITUDE : 121° 4748” W 

MEAN CLIFF TOE ELEVATION (above MLLW): 5.33 M. 
MEAN CLIFF TOP ELEVATION (above MLLW): 8.41 M. 
LONGSHORE OBSERVATION DISTANCE : 1353.00 M. 
BEACH SLOPE : 1:69.56 


SIGMA : 1.26 M. 


TIME INTERVAL RECESSION [X;-X] AVG. EROSION 
[M/YR] a ((M3/YR)/M] 

1940 - 1956 -0.82 0.02 -6.89 

1956 - 1966 -().20 Oo -1.68 

1966 - 1980 -0.63 0.17 -5.29 

1980 - 1984 -3.27 1.93 -20.44 

1940 - 1984 0.84 ©... -7.06 
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Rincon Beach Erosion Results. 
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TABLE VI 
EROSION RATE RESULTS FOR MARINE LAB 


STATION NAME: Marine Lab 

Pa nilUDE 236 47 35° N 

BONGIIUDE = 121° 4732" W 

MEAN CLIFF TOE ELEVATION (above MLLW): 5.15 M. 
MEAN CLIFF TOP ELEVATION (above MLLW) : 8.93 M. 
LONGSHORE OBSERVATION DISTANCE : 506.00 M. 
BEACH SLOPE : 1:66.66 


SIGMA : 1.63 M. 


| [M/YR] oO ((M3/YR)/M] 
1940 - 1966 +1.19 0.45 + 10.63 
1966 - 1970 -2.04 loo -18.22 
1970 - 1976 + 1.30 Oro7 a2 US) 
1976 - 1984 =o 1.10 -12.05 
1940 - 1984 0-4) ee + 4.02 
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Marine Lab Erosion Results (Unsmoothed). 
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Marine Lab Erosion Results (smoothed). 
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EROSION RATE RESULTS FOR NORTH MOSS LANDING 


STATION NAME: North Moss Landing 
LATITUDE : 36° 48’25” N 


LONGITUDE : 121°47’25” W 


MEAN CLIFF TOE ELEVATION (above MLLW): 5.29 M. 


MEAN CLIFF TOP ELEVATION (above MLLW) : 9.81 M. 


LONGSHORE OBSERVATION DISTANCE : 900.00 M. 


BEACH SLOPE : 1:66:66 


SIGMA : 2.08 M. 


TiMetNTeRVAL RECESSION [X;-X] 





[M/YR] G 
1940 - 1970 +0.09 0.03 
1970 - 1976 +2.64 1.26 
1976 - 1980 rOr12 0.07 
1980 - 1984 30 1.56 
1940 - 1984 = 
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North Moss Landing Erosion Results. 


Figure 2.14 
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TABLE VIII 
EROSION RATE RESULTS FOR Zippers DUNES 


STATION NAME : Zmudowski Dunes 

LATITUDE : 36°49’40” N 

LONGITUDE : 121°47’50” W 

MEAN CLIFF TOE ELEVATION (above MLLW) : 3.26 M. 
MEAN CLIFF TOP ELEVATION (above MLLW) : 15.56 M. 
LONGSHORE OBSERVATION DISTANCE: 1050.00 M. 
BEACH SLOPE : 1:69.55 


SIGMA : 0.23 M. 





TIME INTERVAL RECESSION |x;-x] AVG. EROSION 
[M/YR] o ((M?/YR)/M] 

1940 - 1956 -0.16 1.28 -2.48 

1956 - 1970 -0.84 ee -13.07 

1970 - 1984 -0.61 O32 -9.73 

1940 - 1984 052 aa -8.09 
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Zmudowski Dunes Erosion Results (Unsmoothed). 


Figure 2.15 
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Figure 2.16 Zmudowski Dunes Erosion Results (Smoothed). 
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TABLE IX 
EROSION RATE RESULTS FOR MONTEREY BAY ACADEMY 


STATION NAME: Monterey Bay Academy 

LATITUDE : 36° 54'26” N 

LONGITUDE : 121°50’55” W 

MEAN CLIFF TOE ELEVATION (above MLLW) : 4.36 M. 
MEAN CLIFF TOP ELEVATION (above MLLW) : 21.74 M. 
LONGSHORE OBSERVATION DISTANCE : 1220.00 M. 
BEACH SLOPE : 1:76.79 


SIGMA : 0.64 M. 





TIME INTERVAL RECESSION |x;-X| AVG. EROSION 
[M/YR] o ((M?/YR)/M] 
1940 - 1956 a OelS 0.79 ar Ses 
1956 - 1970 =ieo3 1.45 -30.07 
1970 - 1984 -0.03 0.04 -0.60 
1940 - 1984 20033 WI cars. -6.62 
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Figure 2.17 Monterey Bay Academy Erosion Results. 
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lil. WAVES AND TIDES 


A. WAVE CLIMATOLOGY 


An accurate understanding of wave climatology is paramount in any study of 
beach erosion. Wind waves provide the energy to drive the nearshore dynamics 
responsible for changes in coastal morphology. The wave climatology for Monterey 
Bay was calculated using the Wave Information Studies (WIS) deep water variance 
density spectra compiled by the U.S. Army Corps of Engineers. The WIS data is 
broken down into three phases. This study will be concerned with phases | and II 
only. 

In phase I, wind fields are derived from the Gridded Northern Hemisphere 
Pressure Fields of the Navy Fleet Numerical Oceanography Center (FNOC), and the 
North American Historical Weather Map Series of the National Climatic Center 
(NCC). The data from FNOC and NCC represent a twenty year sampling from 1 
January !956 through 31 December 1975. A first guess wind field in regions of low 
streamline curvature is computed from quasi-geostrophic theory. Baroclinic effects 
caused by vertical temperature gradients (thermal wind) are added to the wind 
components and provide the basis for the observed ageostrophic flow. Wind fields for 
synoptic systems exhibiting strong streamline curvature, as typified by an extratropical 
cyclone, are derived from gradient wind theory which accounts for centrifugal 
acceleration in the balance of forces. Additionally, an isallobaric component of the 
wind field is estimated and included in rapidly developing systems where the rate of 
change in pressure with time cannot be ignored. The derived free atmospheric wind 
fields are then reduced to a boundary layer reference level using friction velocity scaling 
(Resio, Vincent, and Corson, 1982). Phase I wind information is computed with a six 
hour time step and interpolated to a 220 KM open ocean grid. In Phase II Surface 
Marine Observations and Airways Surface Observations, both available from NCC, are 
blended into the wind fields computed in phase I, and the modified wind fields are 
interpolated to a 55-km coastal grid. 

Wind waves are hindcast using a discrete spectral wave model developed by 


Resio (1981). Two-dimensional variance density spectra are generated over sixteen 


>I 


directional bands at 0.01l-Hz frequency bands from 0.03 Hz to 0.22 Hz. Dhscrete 
spectral models, in general, are representative for waves in their late stages of 
development; however, they exhibit significant difficulty simulating proper growth in a 
fetch or duration limited domain, or for waves in their early stages of development 
(Dexter, 1974; Resio and Vincent, 1979). Barnet (1968) proposed a more complete 
specification of the classical A + BE(f) source term as a possible solution to this 
problem, where A represents the Phillips (1957) linear, initial wave growth mechanism, 
B , the Miles (1957) exponential, continued wave growth mechanism, and E(f) is the 
energy density spectrum. Although Barnet’s term was three times too low for the 
JONSWAP spectrum, his initiative in incorporating a non-linear input aided in the 
subsequent parameterization of non-linear wave-wave interaction used in the Resio 


(1981) model. The non-linear source term employed by Resio (1981) is given by 
Gay = Da? g* f° expll - (0 1G) (a1) 


From the JONSWAP data, @ = 0.076X, where X is the nondimensional fetch 
X = gx/U? , D is a dimensionless constant, fn 1S the frequency of maximum energy, and 
U is the wind velocity. The linear Phillips mechanism was found to transfer more 
energy into the low frequencies on the spectruin’s leading edge than into the mid-range 
frequencies. It was consequently removed from the Resio (1981) model and replaced 
with the G,) term in the source equation. The total source term for the Resio model 
is then given by the non-linear wave-wave interaction parameterization and the Miles 
exponential growth mechanism 


S = Gi + BE() w(O- 8, ) (3.2) 


ource 


In Equation 3.2, w is a nondimensional directional spreading function where 9, is the 


mean wind direction, and B is given by 
B = 5s f(u cos8)c"! - 0.9] (3.3) 


In Equation 3.3, s is the density ratio between the air and water, and c the phase 


speed, or celerity. 
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Clearly, with the Phillips mechanism missing, another wave generation scheme 
must be included to provide initial wave growth. This initialization. is accomplished 
with a local parametric model. The domain of the parametric model is limited by a 
fixed high frequency cut off which delineates the boundary between the discrete 
spectral and parametric domains. The peak frequency is calculated for each time step 
from equations presented by Hasselmann et. al. (1976). The energy developed in the 
parametric domain is transferred to the discrete spectral domain until the peak 
frequency reaches the high frequency cut-off. The parametric spectrum then shifts to 
an f> equlibrium form, and all growth is continued in the discrete spectral domain. 
The Resio (1981) model provides wave growth and decay rates that are consistent with 
the observed JONSWAP data while maintaining spectral shapes that are highly 
representative of observed spectra (Risio, 1981). In phase I, wave information is 
calculated on a 220-km grid and the spectra are archived at 64 locations in the Pacific 
as depicted in Figure 3.1. In phase II, wave information is calculated using the 
modified wind fields at 55-km grid spacing and the spectra are archived at 53 coastal 
locations as depicted in Figure 3.2. Note that the wind wave model ts a northern 
hemispheric model, and it does not allow wave energy from the southern hemisphere to 


propagate across the model’s boundary. 
B. WAVE REFRACTION 


Wave rays that approach non-parallel isobaths at an oblique angle will tend to 
bend toward the orthogonals of the isobaths as a result of wave refraction. 
Consequently, an area where rays come together is a region of energy convergence, and 
an area where rays separate, a region of energy divergence. In Monterey Bay the Moss 
Landing area at the head of Monterey Canyon provides a good example of energy 
divergence; conversely, Fort Ord typifies a region of ray focusing, or energy 
convergence. 

The effects of shoaling (changing group velocities) and refraction (convergence or 
divergence of variance density) are applied to the WIS deep water spectra to provide 
shallow water and surf zone information. The variance density transformation 1s 
performed using spectral ray theory. In spectral ray theory, one assumes that energy, 
Or Variance, that is associated with a particular frequency and directional band will 
remain in that band throughout the transformation, i.e. constant energy flux 
(Thornton, 1983). Hence the shallow water variance density spectra Sp(F9) are related 


to the deep water variance density spectra S,(f),0)) by 


ais, 


90 


0 
20 


Figure 3.1 WIS Phase ) Grud 
with Archive Locations Indicated. 
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MERCATOR PROJECTION OF PACIFIC OCEAN 


S(68) = K? (68) K? .(£) So(fo,8) J (3.4) 
where Kr is the refraction coefficient, Ks the shoaling coefficient, and J is the Jacobian 
of the transform. 

Here, J = O(f),05)/0(F,8) 
however, since f, = f, J becomes 06(8,)/ (9) and integrates away. 


The shoaling coefficient, Ks, is given from Airy wave theory to first order as the 


ratio of the deep and shallow water group velocities. 


B= 6 (Cpieal (3.5) 


where C g is given by 


Cy = C/2(1 + 2kh/sinh 2kh) (3.6) 
In shallow water, surface gravity waves become non-dispersive and the group velocity 


can be approximated by 
Cy Se (Gy) (3.7) 


The deep water group velocity is fixed for a particular frequency, whereas the shallow 
water group velocity will change with shoaling water depths; hence Ks will increase as 
the wave crest approaches the shore. 


The refraction coefficient, Kr, provides a measure of ray separation 
Kr = a/(b,/b) (3.8) 


where by is the deep water ray separation, and 6 is the shallow water ray separation. 
Values of Kr that are greater than unity represent energy convergence, and those less 
than unity, represent energy divergence. The refraction coefficients for this study were 
calculated using the linear refraction model by Dobson (1967) modified to refract rays 


from a specified shallow water location, back to deep water. This “back” refraction 
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model was developed to obtain refractive information at specific locations. The model 
searches for the 4-m isobath along a normal to the shoreline that originates at an 
observation site or other point of interest. Rays for each frequency in the WIS 
spectrum (0.03 Hz - 0.18Hz) are started at a specified location and depth, nominally 
4-m, and are propagated offshore at 0.l-degree increments over the range of all 
possible arrival angles. Once in deep water, the rays are stopped and the deep water 
angle is measured; then they are turned around and propagated back inshore along the 
same path to calculate the spectral refraction coefficient. The rays must return to the 
initial location within the specified area of £100 m to be a valid calculation. 
Examples of spectral refraction in Monterey Bay are presented in Figures 3.3, and 3.4. 
The bathymetry used in the refraction model is critical to the analysis. The 
accuracy of the calculated refraction information can be no better than the accuracy of 
the bathymetry; consequently, considerable effort was devoted to accurately depict the 
bottom. Original NOAA data was obtained that had been projected onto a six-second 
modified Universal Transverse Mercator (UTM) grid. The bathymetry data was 
initially screened for bad points, then the data was scanned along meridians and 
parallels for changes in slope that exceeded 30°. Points that generated unrealistic 
bathymetry were extracted from the data base, and the resultant bathymetry was 
projected -onto an x-y plane via a modified UTM projection. The bathmetry was 
interpolated to a 200-m rectangular grid using a piecewise-linear, triangular plane 
interpolator (Franke, 1986). The triangulation method provided reasonable results in 
data sparce regions, and minimum distortion of bathymetric features. Intermediate 
smoothing was accomplished using a nine point weighted linear averager, and further 
smoothing was provided by the model itself which calculates bottom curvature by 
fitting a quadratic surface to adjacent isobaths. Waves originating from the northerly 
most quadrant are refracted when travelling over the shoal shelf region between the 
Farallon Islands, and Point Santa Cruz. To include the refraction to the north, a 
supplemental Northern bathymetry was gridded for a section of the California coast 
north of Santa Cruz. Refraction coefficients were calculated for the portion of the ray 
path which traversed the northern bathymetry, and for the portion of the ray path in 
Monterey Bay. The two refraction coefficients were multiplied together to determine a 
complete refraction coefficient for use in the spectral transformation. The Monterey 


Bay bathymetry is shown in Figure 3.5. 


a7 





Figure 3.3. Ray Trace for Marina, f = 0.05 Hz.. 
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The objective in performing the spectral transformation is to produce a wave 
height that can be used to calculate wave set-up, run-up, and energy. The significant 


wave height, Hs, is related to the total energy, or variance, o* , by 
ety Adena 4 2V (8G) “ae (3.9) 


The variance is calculated by integrating the transformed shallow water spectrum over 


all frequencies and directions. 


o* = fg Sy, (£8) df dé (3.10) 


C. WAVE SEY-UP AND RUN-UP 


As waves break in the surf zone, they transfer momentum shoreward. The 
change in shoreward momentum flux is primarily balanced by hydrostatic pressure; 
accordingly, the hydrostatic pressure balance is maintained through an increase in the 
mean sea level shoreward of the surf zone. This balance assumes that there are no 
mean currents, and no longshore gradients in the shoreward momentum flux, and 1s 


given by 
6/Ax{2EC,/C - E/C] = -pg(n +. hjén /éx | (Saials) 


where the term in brackets on the left hand side is the momentum flux due to the 
normally incident breaking waves. As waves break, the wave energy, E, decreases 
resulting in a decrease in momentum flux, 1.e. negative gradient, which must be 
balanced by a positively sloping sea surface. This mechanism is known as set-up (n) 
and its maximum height above mean sea level occurs at its shoreward limit. Wave 
set-up studies by Guza and Thornton (1981) suggest that wave set-up 1s proportional 
to the deep water RMS wave height. Additional work by Holman and Sallenger (1985) 
over a wider range of wave and beach parameters, corroborate the results of Guza and 
Thornton, and further suggest that there is a correlation between the tidal stage, and 
wave Set-up. Beach erosion primarily occurs at high tide, and consequently this study 


uses wave set-up at high tide (Holman and Sallenger, 1985). 


5g 





Figure 3.4 Ray Trace for M.B. Academy, f = 0.09 Hz.. 
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1) = 0.88H’ (3.12) 


The significant wave height, H, is calculated at the 4-m isobath, where H Pee O eect 
the wave is back shoaled without refraction to obtain a deep water significant wave 


height, H’., , using 
Hg = Hs/Ks (3.13) 


where Ks is calculated using the peak frequency. Wave set-up is dipicted in Figure 3.7. 

Wave run-up is the maximum superelevation of the sea level caused by the rush 
of breaking waves up the beach. Thornton and Guza (1981) suggest that run-up is also 
proportional to the deep water wave height H'.g - Holman and Sallenger (1985) note a 
tidal correlation and find run-up at high tide to be given by 


nN = 0.14H’ (3.14) 


Wave set-up and run-up are important parameters in raising the mean sea level to a 


stage where waves can erode the cliff base. 
D. TIDES 


The calculation of tidal height is inherently critical to the specification of mean 
sea level. The simultaneous presence of high tide and storm waves presents the highest 
probability that an erosion event will occur. Consequently, any parameterization of 
the fluxuating sea surface must include a harmonic tidal analysis. The harmonic 
analysis allows one to express the tidal height at any time, t, as a finite sum of its 
harmonic constituents, each of which is representative of a specific astronomic 


condition. The tidal expression from Dennis and Long (1971) 1s presented below 
ae ol as ee Pid COsmd sec vga l,)| Gyiley) 


where Hp is the height of the mean tidal level (above datum), f,, is the node factor of 
the constituent, H, is the mean amplitude of the constituent, a, is the speed of the 
constituent, t is the time reckoned from some arbitrary point, KS is the epoch of the 


constituent at t = 0 for the observation period, [V) + uJ, is the value for the 
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Figure 3.6 Bathymetry of Monterey Bay. 
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argument at t = O for the observation period, and n is the particular constituent being 
computed. 

The constituents for Monterey Bay were obtained from a harmonic tidal analysis 
performed by the Tidal Predictions Branch of the National Ocean Service. Their data 
base was a thirteen-year time series of tidal data observed at Monterey’s Municipal 
Wharf No. 2. Twenty constituents were isolated with the threshold criterion that a 
constituent would be consideried valid only if its amplitude contribution was greater 
than 0.6 cm. The speed, node factor, and argument of the constituents were extracted 
from Schureman (1971). The twenty constituents were used to calculate tidal heights 
at every hour during the twenty-year WIS time series. 

The tides in Monterey Bay are typically semi-diurnal with two low and two high 


water cycles per day. The mean tidal range is 0.98 m. 
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Figure 3.7 Wave Set-up. 
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IV. AN EROSION MODEL FOR MONTEREY BAY 


A. MODEL PRESENTATION 


The objective of this thesis is to develop an erosion model for Monterey Bay. 
Lima and Sklavidis (1984) computed a multiplicative constant, k, for the Brunn erosion 


model presented in Komar (1983) 
R = K((S - T)/tanB] (4.1) 


where R is the recession rate, S is the height of mean sea level, T is the cliff toe 
elevation, tanB is the beach slope. The value k = 0.000098 was calculated from a 
linear regression of the Lima and Sklavidis (1984) recession rate data to the Brunn 
model. — 

This research presents a model where recession is a non-linear function of the 
deep water significant wave height H’,, (for purposes of clarity in the following 
equations, H’., will be replaced by H). 


Rmod = [AH* + BH + C(S - T)//tanB (4.2) 
Rad is the modeled recession rate and A, B, and C are constants. The deep water 
significant wave height is determined as discussed in Chapter III and the mean sea 


level, S, is given by 
S = Tide + Set-up + Run-up (4.3) 


The beach slope, tanB, was extracted from large scale hydrographic charts and varied 
little along the Bay. The idea of the model is that bluff recession only occurs when 
superelevation of the mean sea level, caused by tides and waves, exceeds the height of 
the cliff toe. When this criterion has been met, recession will occur at a rate which is 
proportional to: the wave energy, All?, the significant wave height, BH, and the 


ammount of water which exceeds the cliff toe, C(S-T). 
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In solving for A, B, and C, it is necessary to minimize 


¥ (Rmodg > Robs)” = MIN (4.4) 


where R are the recession rates observed during the time-span of the WIS spectral 


obs 
data from 1956 to 1975. The partial derivatives of equation 4.4 are taken with respect 


to A, B, and C, and set equal to 0 to define the following three normal equations 


ADH) BH? + OR(S-TIA* = Tani Ree (4.5) 
AYH?® + BYH? + C)(S-T)H = TanB Y}Ro,.H (4.6) 
AYH*(S-T) + BYH(S-T) + C)\(S-T)* = TanB YRop,(S-T) (4.7) 


The three equations are then solved simultaneously by the matrix equation 


yak BE Y(s-DH2| [A) | Tanp TRopsH? (4.8) 
YH? YH? Y(S-T)H |X | B | = | TanB YRop.H 
YH*(S-T) YH(S-T) }¥(S-T) C}  |TanB }’Rop,(S-T) 


This procedure is repeated for each of the twelve stations and the resultant coefficients 
are averaged to determine a mean A, B, and C for Monterey Bay. 

The model is very sensitive to changes in the elevation of the cliff toe, as this 
elevation is the model’s go, no-go threshold value. An initial mean sea level model was 
run at each station to determine the number of hours, N, that the sea height exceeded 
the cliff toe for the twenty year WIS data. Adjustments , less than 0.5 m, were then 
made to the surveyed cliff toe elevation to ensure that the modeled N for each station 
was physically reasonable, and its value was within one order of magnitude of the other 


stations. The adjustments themselves are not unreasonable considering the inherent 
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variability of the cliff toe elevetion due to falling cliff material, and seasonal sediment 


transport. 


B. DISCUSSION OF RESULTS 


Recession model coefficients A, B, and C are determined for all sites that exhibit 
recession between 1956 and 1975. Moss Landing Marine Laboratory, and North Moss 
Landing are the only sites where continuous accretion during the time period prevented 
application of the model. Initial variance-and standard deviation values are determined 
for the coefficients at the remaining ten sites, and based on this initial analysis, the 
coefficients for Monterey Sand Co. are rejected as being statistically unrepresentative 
of the erosion processes in Monterey Bay. The coefficients of the remaining nine sites 
are averaged to provide a mean A, B, and C, and hence define the recession model. 
Wiese coetiiciemts are listed in Table X. 

The calculated mean values of the coefficients provide a reasonable depiction of 
bluff recession for regions of low wave energy, and will heretofore be referred to as the 
Low Energy Model. The linear wave height coefficent B, is the most significant, and 
defines a mildly sloping recession profile. Unfortunately, as can be seen in Figure 4.1, 
the Low Energy Model becomes unreliable in regions of energy convergence. A 
significant wave height of 5 m appears to be the threshold where the model’s energy 
term, AH? , begins to have a significant impact on the model recession. Hence a new 
value, Asp , for the energy coefficient, A must be determined from a weighted average 
of the A coefficients computed in the following regions of high wave energy: South 
Fort Ord, Stillwell Hall, Marina, and Monterey Bay Academy (see Table X). A FEligh 
Energy Model is then obtained by replacing A with Ajyyp . The coefficients B and C 
remain the same as in the Low Energy presentation, however, the new, Ayjp 
parameterization of the modelled energy produces a much more realistic model 
response in energetic regions. This improved response is also depicted in Figure 4.1. 
The high energy recession curve is too shallow at low wave heights and consequently, 


the model overpredicts the observed recession in these regions. 
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Station Name 


NPS Beach Lab 


Phillips Petroleum 


Sand City 


Monterey Sand Co. 


South Porm Ord 


Stillwell Hall * 


Marina * 


Rincon Beach 


Zmudowski Dunes 


Monterey Academy * 


Mean Values 


High Energy values * 


Standard Deviations 


* Indicates regions of high wave energy, He ESS OVA G0 


TABLE X | 
CALCULATED RECESSION COEFFICIENTS 


A [1/M?] 


-0.058619 


-0.025892 


-0.052943 


-0.063531 


-0.011144 


-0.023181 


-0.016920 


-0.010737 


-0.020776 


-0.023647 


-0.027095 
-0.02005 1 


= 0.0287 
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B [1/M] 


0.233661 


0.14414] 


0.358974 


0.461977 


OxvlS3575 


O'235529 


0.196066 


0.101322 


0.188279 


0.228347 


0.20203 


same 


£ 0.214 


C {1/M] 


-1.79972 x 10° 
-2.29962 x 10° 
-2.59265 x 10°9 
1683579 Oue 
5.18 795eL0m 
2.43953 x 1077 
-1. 205480. ule 
5.60665 x 10> 
-2.08930 x 1073 
3.02067 x 10°? 


-0.002317 


same 


+ 0.0024 


RECESSION MODEL 


SIG 


\ 


8.0 


V4 
BP 


IN 


A 


ip 


e 
. 
‘< 


© =) © © em 
IN (oO Oe. a: BD N 


Coca re mee anon 1S 


Figure 4.1 Recession Model Profile. 
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C. CONCLUSIONS 


The recession rates are found to be related to the local wave height. A 
polynomial regression model in wave height shows the recession to be most correlated 
with wave height to the first power (B is the largest coefficient) and negatively 
correlated with wave energy, H*. The recession observations at Sand City and 
Monterey Sand Company show significantly more bluff recession than the Low Energy 
Model suggests, and the observed significant wave height does not warrent use of the 
High Energy Model. Local sand mining operations are the most probable cause for 
the high recession values at these sites. For the remaining eight stations, however, the 
most representative results are obtained by applying the Low Energy Model in regions 
of low wave energy, and the High Energy Model in regions of high wave energy. In 
doing so, the standard error calculated between the observed and the computed 
twenty-year recession rates is £3.41 m which equates to +0.17 m/yr. The mean 
recession rate for the stations which experienced recession between 1956 and 1975 , 
discounting stations Sand City and Monterey Sand Company, 1s 0.76 m/yr, hence the 
model represents actual recession to +22 %. These results are listed in Table XI. 
Graphic depictions of the twenty year recession at each site are depicted in Appendix 
C. The observed recession results presented in Table XI only represent events occuring 
between January 1956 and December 1975, and a complete summary of bluff recession 
between 1940 and 1984 is provided in Table XII. 
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TABLE XI 


RECHESSIONMMODEL RESULTS (1956-1975) 


Station Name 


NPS Beach Lab 
Phillips Petroleum 
Sand City 
Monterey Sand Co. 
South Fort Ord 
Stilwell Hail 
Marina 

Rincon Beach 
Zmudowski Dunes 


Monterey Bay Academy 


The Standard Deviation for the Model is + 3.41 m for 20 


Observed 


9.00 m 


7.80 m 


23.00 Ti 


Sse) 00 


Heieey nal 


19.40 m 


22.4 m 


9.60 m 


16.80 m 


21.8 m 


years which equates to + 0.17 m/jr. 


i 


Modelled 


11.45 m 


13.43 m 


14.4 m 


14.86 m 


10.57 m 


10.97 m 


e350 Tl 


io ml 


14.08 m 


iS. mM 


High Energy Model 


19.68 m 


20.49 m 


23.49 m 


19.97 m 


TABIEE xh 
OBSERVED RECESSION SUMMARY 


Station Name Recession (m/yr) 
NPS Beach Lab -0.58 
Phillips Petroleum -0.85 
Sand City -0.97 
Monterey Sand Co. -1.94 
South Fort Ord -1.12 
Stillwell Hall -1.98 
Marina -0.38 
Rincon Beach -0.84 
Marine Lab + 0.45 
North Moss Landing + 0.03 
Zmudowski Dunes -0.52 
Monterey Bay Academy -0.33 
Standard Deviation +0.73 Mean = -0.75 m/yr 
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APPENDIX A 
PROGRANI LISTING PHOTODAT 


= INPUT “WHEAT IS yOuU2 STSTION 2.D. 1.8. LABISSE":Es 
3 GEEN £3 FOR DUTELT AS a: 
4 2RINT 
Epoeiit. = bas 
7 PRINT 
8 COUNT=2 
9 PRINT “REMEMEE2 79D ALWAYS AOE RCACH YOUR FOINT FROM THE SAME DIRECTION" 
12 PRINT "TO PREVENT SLESTRONIC OR MECHANICAL BACKLASH FROM OCCURING IN" 
11 PRINT "THE DATA. THE FIRST STED IS TO DEFINE YOUR ORIGIN. MOVE THe" 
12 PRINT "VIEWER SVER TRE CORRECT POINT AND PRESS (¢ DH ) (RPETURN) AT TRE" 
1S PRINT "(€ : ) PROMPT. REDEAT THIS PROCEDURE FOR EACH AXIS." 
14 ERINT 
1S sengee2 
16) SeNa 
fe 2e0N7 
19 PRINT “REDEFINE THE CRIGIN, PRESS ( Dn) AT THE FROMET" 
ei F@R I=1 TO 3 Ssre> 2 > THIS LOGE INITIALIZES TKE X&Y AXES 
SS J=I-1 1 N= THE DMC-400 READ/WRITE ADDRESS 
32 N=100B+JI 7 N=100@0 X-AXIS N=19@2 Y-AXIS 
3S as: > AXIS COUNTER. X=1 Y=2 
41 PRINT CHRSCINEON)): 
45 B=IND(N+1) > CHECKS DMC-4a2a STATUS ADDRESS 
SQ@ C=HINT(s/2 7 "CY" CHECKS STATUS FLAG HIT 
Sami= G=S/72 Tees 62 > DOES DMC-4@a@ 4AVE DATA TO SEND? 
6@ Ie INT(C/2) =C7=E THEN 45 ’ IS DmC-40a READY TO ACCEPT. COMMANDS? 
65 2= a=4 THEN 3S 
72 AS=INKEYS 1 INPUT COMMAND FROM KEYBOARD 
7i IF LEN(AS) =@ TREN 45 
7=z PRINT As: 
SOT Ny ASCKASs) 1 SEND THE “DEFINE HOMES" COMMAND 
7 AsxsG-j 
78 GET 4S 
BQ AS=THRSE CINE IND | + READ DATA FROM THE Dmc-4au 
S IF GBi() CRRS(22) TREN ORINT AB? 
92 IF GAs=CHRS(65, “HEN 18 
94 GOTO 45 
95 NEXT : 
1a PRINT "LEAVES Tee «/Y CONTROLS ALONE AND CHECK TO INSURE THAT x & Y =a" 
lige 2RINT 
195 INPUT "SRESS ¢ & ) WHEN READY" :Fs 
L112 IF F$s"0" TREN 149 
1:5 GOTS 105 
129 ORINT 
LS2 INGWT “HAS "R= ORIGIN BEEN CORPECTLY DEFINED? (Y/N) "29S 
L335 IF Dss"Y¥" TREN 145 
242 GOTO :6 
144 PRINT 
145 PRINT “YOU P= NGw READY TO BEGIN TAKING DATA. MOVE THE VIEWER OVER" 
146 SRINT "TRE DTESI2ED AOSITION AND ORESS ¢( DH ) TO CATALOG TKE FCOINT. START" 
147 ORINT "WIT Tea BASELINE, THEN SCALING POINTS. THEN THE SHORELINE ZOINTS. " 
148 CRIN 
i199 I=1 
eviad =OR I=1 TO S STEP = * LOOE TG TAKE xX/7Y DATE 
213 J=I-1 
ESQ V=1 92+ 5 
SSz DRINT 
225 3SPITK™ "S9SSS ¢ ~> ) “OG CS7-ALCS "Ni" ADDRESS DATA, " 
22a GCSLe Sev 7 SUBROUTINE TO COMMUNICATE w/DMC-saa 
a — a, rd "Aye ToS toa 
ea aa 
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OOO YNNNN NEN IND 
Swoon gwar ola 


fr 0) fo a Fa fu a a TO a A fat Pe fs 
Gy fia ee 


) ee ad 

Se.0O 2748 

SF ON=1202 "HEN 272 

GCTO 75 

Y=Re (-1) 

GOTO 278 

Z=& 

NEXT I 

SRINT 

PRINT "YOUR x, Y & Z COORD’S ARE ";X,Y,Z 

INPUT “WOULD YOU LIKE TO STORE THESE PTS IN THE DATA FILE (Y/N) :":Y% 
IF Y$="N" THEN 289 

PRINT #1,X:¥32Z 

COUNT=COUNT+1 

IF X=@ THEN 129 

IF Y=@ THEN 129 

PRINT 

INPUT "WOULD YOU LIKE TO CATALOG MORE DATA POINTS? (Y/N) ":G6$% 

IF GS0"N" THEN 148 

PRINT 

PRINT “THERE ARE ";sCOUNT:" POINTS STORED IN FILE ";Es 

PRINT 

PRINT “BE SURE TO LOG THIS INFORMATION, AS YOU Witmeneaern tyemeneR «* 
PRINT 

PRINT "THAT’S ALL FOLHS!'! ITS MILLER TIME !!'" 

CLOSE 

END 

REM THIS IS A SUBROUTINE TO COMMUNICATE BETWEEN THE DMC-42a8 5s IBM PC 
A=1 1 THIS SUBROUTINE COMMANDS THE Dmc-4aa 
m=6 . TO REPORT ITS FOSITION (TELL FOSITION 
Q=2 7 COMMAND), AND THEN READS THE POSITION 
x= 1 7 FROM THE N ADDRESS. THE POSITION (CIN 
PRINT CHRSC(INP(N));: 

B= INP (N+1) 

C=INT(S/2) 

IF C=B/2 THEN S65 1 CHECKING STATUS & RERDING DATA 
IF INT(C/2) =C/2 THEN S22 

IF Q=4 THEN 692 7 AXIS COUNTER 1=X =Y 
AS=INKEYS 

IF LEN(AS)=0 THEN S20 7 ASSIGN "TELL POSITION COMMAND" 
PRINT AS; 

OUT N, ASC (AS) ’ SEND COMMAND TO DmMC-4aa 

A=A+1 

GOTO Sza 


AS=CHRS( INP CN) ) 
TE eae) GaSe) 
IF AS=CHRS (ES) 


TREN PRINT AS: 
Tae Neo we 


‘VeASC (As) 

IF (V)=48) AND (V¢é=57) THEN 6e2 ’ 
IF (V)=65) AND (V¢«=720) THEN G22 ’ 
GOJOgss 


IF V¢6éil THEN Z=(V—-48) ? 
IF V¥) 60 THEN Z=<V—sSS5 

R=R+ (Z*1E6% (M—1)) 

IF Msé THEM Ze=Z 


m=M—1 

Me = alas) Ses") 
K=K+1 

SOveeas 

te ce. < 8) THEN D7 1G) 


R=R-167772E16# 
RETURN 
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CHECAING JGese= l- Danas ORe 
CHARACTERS 8-9 OR A-F 


CHANGING FROM REX 10 DEE 


12 
2B 
52 
40 


i 
-_ 


sa 

72 

82 

92 

98 

1a 
Las 
Lia 
ives 
152 
162 
172 
198 
2a 
F212 
ee2 
23 
=42 
245 
=48 
E252 
e72 
z74 
AP ga 
ae ial 
= 8 
oc 
265 
eo 
29a 
272 
eJe 
S22 


rand 
we 


322 
23a 
Sao 
362 
4&2 
420 
452 
462 
S22 
6B 
Wa 


APPENDIX B 
PROGRAM LISTING WRKPROG 


INFLT "“INGUT THE NUMBER OF POINTS IN THE DATA FILE:":N 

DIM X#(N).Y#(N), Z#(N), XR#(N), YR#ON) 

PRINT 

INPUT “INEUT THE NAME OF YOUR DATA FILE. I.E. LABISS3S:"3ES 
PRINT 

OPEN “I",#1,89 

[oR l=. 70 N 

INPUT #1,X#(I),Y#C(I), Z#(I) 


NEXT I 
I=1 
AX#=, AWJVL2S762875# . REGRESSION COEFFICIENTS TO CONVERT 
BX#=-. BOWB7EZ47924A4H , ENCODER COUNTS TO CENTIMETERS. WHERE 
AY#=. ANABGISESSIL7S# ee =a Coe owing) -. 
BY#=-, 095152291474% 
FOR T=2 TON 
X# (IT) =X8#(1I) #QX# + BX# * CONVERSION E-COUNTS TO CM 
Y#(I)=Y#(I) #QY# + BY# 
NEXT I 
ALFA#=ATN (Y# (2) /X#(2) ) * COORDINATE SYS ROTATION ANGLE 
FOR J=1 TON 
XR#¥( J) =(X#(5) COS (ALFA#) )+(VY#H(J) #SIN(ALFAS#) ) * COORDINATE TRANSFORMATION 
YR#(5) = (0 (-1) *X8(J) #SIN(CALFA#) ) + (Y#( SI) #COS (ALFA) ) 
NEXT J 
PRINT 
INPUT “INPUT THE MEASURED FIELD DISTANCE CIN METERS) :"s:GDIST 
PRINT ; 
PDIST#=SOR ((XR#(4) -XR#(3))72 + (YR (4) -YR#B(S) 72) 
or Ny 
LFERINT “SHOTODISTANCE = "“sPDIST#s;" Cm" 
PRINT 


PSCALE#= (FP DISTH#/1A0) /GDIST 
PRINT FPSCALE# 


Pe=1/SSCALE# 
fea tt 
LERINT “SHOTOSCALE IS: “;P# 
LPRINT 


LPRINT “YOUR ACTUAL SCALE, ROTATED X/Y PAIRS FOLLOW (METERS) " 
POS H= 3 TO N 

XR#CK)P = (XR#(K) S/PSCALE#) /122 

YR#(K)= (YR#(K) /PSCALE#) /100 

LERINT XR#(K). YR#(K) 

Nex Kh 

LERINT 

pogece= 2 10. (N—1) * TRAPAZOIDAL RULE INTEGRATION 
AREA#=AREA# + ((CYR#(L) +YR#(L+1))/2) #ABS (XR#(L+1)—-XR#(L) ) 

NEXT L 

PRINT 


LERINT "THE AREA IS: "“s;AREQA#;" SQUARE METERS FOR STATION ";ES 
CLOSE 
END 
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APPENDIX C 
RECESSION MODEL RESULTS 
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